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Abstract. We consider the Navier-Stokes system in a bounded domain with a smooth boundary. 
Given a sufficiently regular time-dependent global solution, we construct a finite-dimensional feedback 
control that is supported by a given open set and stabilizes the linearized equation. The proof of 
this fact is based on a truncated observability inequality, the regularizing property for the linearized 
equation, and some standard techniques of the optimal control theory. We then show that the control 
constructed for the linear problem stabilizes locally also the full Navier— Stokes system. 
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1. Introduction. Let C M"^ be a connected bounded domain located locally 
on one side of its smooth boundary F = d^l. We consider the controlled Navier-Stokes 
system in fl: 



Here u ~ (ui, U2, U3) and p are unknown velocity field and pressure of the fluid, i> > 
is the viscosity, {u ■ V) stands for the differential operator uidi + U2d2 + us^s, h is 
a fixed function, and C is a control taking values in the space £ of square-integrable 
functions in J7 whose support in x is contained in a given open subset C 17. The 
problem of exact controllability for (|1.2I) was in the focus of attention of many 

researchers starting from the early nineties, and it is now rather well understood. 
Namely, it was proved that, given a time T > and a smooth solution u of 
(|1.2|) with C = 0, for any initial function ug sufficiently close to u(0) one can find a 
control C, : [0, T] — >■ 5 such that the solution of problem (|1.2|) supplemented with 

the initial condition 



is defined on [0,T] and satisfies the relation u{T) = u{T). We refer the reader to [6l 
[TT| [T2I [T3l [T4| [7] for the exact statements and the proofs of these results. 

Even though the property of exact controllability is quite satisfactory from the 
mathematical point of view, many problems arising in applications require that the 
control in question be feedback, because closed-loop controls are usually more stable 
under perturbations (e.g., see the introduction to Part 3 in [5]). This question has 
found a positive answer in the context of stabilization theory. It was intensively 
studied for the case in which the target solution u is stationary (in particular, the 
external force is independent of time). A typical result in such a situation claims that, 
given a smooth stationary state ii of the Navier-Stokes system and a constant A > 0, 
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dtU+{u-\/)u-i^Au + \/p = h + C, V-it = 0, 
u|p = 0. 



(1.1) 
(1.2) 



u(0, x) — uq(x) 



(1.3) 



1 



2 



V.BARBU, S.S.RODRIGUES, AND A.SHIRIKYAN 



one can construct a continuous linear operator : — > £ with finite-dimensional 
range such that the solution of problem - (I1.3P with C, = Ku{u — u) and a 

function uq sufficiently close to w is defined for all t > and converges to u at least 
with the rate e~'^*. We refer the reader to the papers [9l [lOl |T7l [181 [TH [1] for boundary 
stabilization and to [U |4l [3] for stabilization by a distributed control. 

The aim of this paper is to establish a similar result in the case when the target 
solution u depends on time. Namely, we will prove the following theorem, whose exact 
formulation is given in Section HI 

Main Theorem. Let {u,p) be a global smooth solution for problem (jl.ip . (|1.2p 
with C = such that 

ess sup \&ld"u{t,x)\ < R for j = 0, 1, |a| < 1, 
{t,x)GQ 

where Q = R+ x f2 and R > is a constant. Then for any A > and any open 
subset UJ d there is an integer M = AI{R, X,uj) > 1, an M-dimensional space 
£ C C(j"(a;, M'^), and a family of continuous linear operators Ku{t) : L^(fi,K."^) — > E, 
t > 0, such that the following assertions hold. 

(a) The function t i— >■ Kji(t) is continuous in the weak operator topology, and its 
operator norm is bounded by a constant depending only on R, A, and lo. 

(b) For any divergence free function uq £ HQ{fl,]S?) that is sufficiently close 
to u{0) in the H^-norm problem (jl.ip - (|1.3p with ( = Ku{t)(u — u(t)) has a 
unique global strong solution {u,p), which satisfies the inequality 

\u{t) ~ u{t)\H^ < Ce-^'\uo - u{Q)\m , t>0. 

Note that this theorem remains true for the two-dimensional Navier-Stokes sys- 
tem, and in this case, it suffices to assume that the initial function uq is close to u(0) in 
the L^-norm. Furthermore, the approach developed in this paper applies equally well 
to the case when the control acts via the boundary. This situation will be addressed 
in a subsequent publication. 

As was mentioned above, the problem of feedback stabilization is rather well un- 
derstood for stationary reference solutions. Let us explain informally the additional 
difficulties arising in the non-stationary case and reveal a common mechanism of stabi- 
lization. A well-known argument based on the contraction mapping principle enables 
one to prove that a control stabilizing the linearized problem locally stabilizes also 
the nonlinear equation. Thus, it suffices to study the linearized problem. The main 
idea in the stationary case is to split it into a system of two autonomous equations, 
the first of which is finite-dimensional and has the zero solution as a possibly unstable 
equilibrium point, whereas the second is exponentially stable due to the large negative 
eigenvalues of the Laplacian. One then applies methods of finite-dimensional theory 
(e.g., the pole assignment theorem [21] Section 2.5]) to find a stabilizing feedback 
control for the first equation and proves that using the same control in the original 
problem yields an exponentially decaying solution; see [2J HI [T71 131 HB HH [I] . 

It is difficult to apply this approach in the case of time-dependent reference so- 
lutions, because a non-autonomous equation does not necessarily admits invariant 
subspaces. However, the above-mentioned scheme for stabilization is based essen- 
tially on the so-called Foia§-Prodi property for parabolic PDE's [8]. It says, roughly 
speaking, that if the projections of two solutions to the unstable modes converge to 
each other as time goes to infinity, then the difference between these solutions goes 
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to zero. It turns out that the conclusion remains true if the projections arc close to 
each other at times proportional to a fixed constant. The main idea of this paper is to 
choose a control that ensures the equality at integer times for the projections of two 
solutions to the unstable modes. More precisely, we consider the following problem 
obtained by linearizing (jl.ip . (|1.2p around a non-stationary solution u(t,x): 

dtv + {u-V)v + {v-V)u-iyAv + Vp = C, V-i; = 0, u|r = 0. (1.4) 

Let us assume that, for a sufficiently large integer N, we have constructed a continuous 
linear operator ^ : L^{il, M"^) — > i^((0, 1); £) such that, for any initial function vq, the 
solution of (|1.4p with C ~ ({"^o) issued from vo satisfies the relation 11^^(1) = 0, 
where IIjv stands for the orthogonal projection in onto the subspace spanned by 
the first N eigenfunctions of the Stokes operator in fl. In this case, using the Poincare 
inequality and regularizing property of the resolving operator for (|1.4p , we get 

\v{1)\l2 = \{I- Un)v{1)\l^. < Cia~'/'K.(l)|Hi 

< C2a^^^^(|vo|L2 + |C(wo)|l2((o4);£)) < Csa^^^^lwoU^, (1.5) 

where {oij} denotes the increasing sequence of the eigenvalues for the Stokes operator 
and d, i = 1,2, 3, are some constants not depending on N. The fact that C3 is inde- 
pendent of iV is a crucial property, and its proof is based on a truncated observability 
inequality (see Section [5731 in Appendix). It follows from (jl.5p that, if is sufficiently 
large, then |u(l)|i2 < e~'^|uo|L2- Iterating this procedure, we get an exponentially 
decaying solution. Once an exponential stabilization of the linearized problem (|1.4p is 
obtained, the existence of an exponentially stabilizing feedback control can be proved 
with the help of the dynamic programming principle. We refer the reader to Section [3] 
for an accurate presentation of the results on the linearized equation and some further 
comments on the existence of Lyapunov function, derivation of a Riccati equation for 
the feedback control operator, and the dimension of controllers. 

The paper is organized as follows. In Section[2l we introduce the functional spaces 
arising in the theory of the Navier-Stokes equations and recall some well-known facts. 
Section |3] is devoted to studying the linearized problem. In Section 21 we establish the 
main result of the paper on local exponential stabilization of the full Navier-Stokes 
system. The Appendix gathers some auxiliary results used in the main text. 

Notation. We write N and M for the sets of non-negative integers and real num- 
bers, respectively, and we define No = N \ {0}, M.^ = (0, +00), and M^. = (s,-|-oo). 
We denote by C M"^ a bounded domain with a C^-smooth boundary F ~ dil, and 
for r e R, we set = {t,t + 1), Q^. = 1^. x fi, and Q = R+ x O. The partial time 
derivative dtu of a function u{t, x) will be denoted by Uf. 

For a Banach space X, we denote by | • the corresponding norm, by X' its 
dual, and by (•, ■)x',x the duality between X' and X. 

If X and Y are Banach spaces and / C M is an open interval, then we write 

W{I, X, Y) := {/ e L\I, X)\ fte L\I, Y)}, 

where the derivative /t = is taken in the sense of distributions. This space is 
endowed with the natural norm 

\f\w(i,x,Y) {\f\\^{i,x) + l/tli2(7.y)) ^ ■ 
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Note that if X ~Y, then we obtain the Sobolev space W-'-'^{I,X). 

If / C M is a closed interval, then C{I,X) stands for the space of continuous 
functions f : I ^ X with the norm 

\f\c{i,x) =max|/(t)|x. 

For a given space Z of functions / — f{t) defined on an interval of M and a 
constant A > 0, we define 

Zx := {.feZ \ e(^/2)V € Z}. 
This space is endowed with the norm 

Throughout the paper, we deal with two integers, TV and M. Roughly speaking, 
N stands for the number of unstable modes in the linearized Navier-Stokes system 
and M denotes the space dimension of the control function arising in various problems. 

C'[ai,...,afc] denotes a function of non-negative variables Uj that increases in each 
of its arguments. 

Ci, i = 1, 2, . . . , stand for unessential positive constants. 

2. Preliminaries. 

2.1. Functional spaces and reduction to an evolution equation. In what 
follows, we will confine ourselves to the 3D case, although all the results remain valid 
for the 2D Navier-Stokes equations. 

Let r2 C M'^ be a connected bounded domain located locally on one side of its C^- 
smooth boundary F = dfl. It is natural to study the incompressible Navier-Stokes 
system as an evolution equation in the subspace H of divergence free vector fields 
tangent to the boundary: 

H — {ue L'^{n,R^) I V - u = in fl, u - n = on F}. 

Here i^(ri, K^) is the space of square integrable vector fields (ui, U2, u^) in Jl, V • w := 
+921*2 is the divergence of u, and n is the normal vector to the boundary F. 
Let us denote by H''{il) the Sobolev space of order s and by iJ'*(f2,K^) the space of 
vector fields in J7 whose components belong to H''{n). To simplify notation, we will 
often write and i?**; the context will imply the domain on which these spaces are 
considered. Define 

V ■.= {ue H\n,R^) I V • u = in 17, u = on F}, [/ := H^{n,R^)nV. 

Note that U coincides with the natural domain D{L) of the Stokes operator L = 
—vHA, where 11 is the orthogonal projection in L^(ri,R'^) onto H. The spaces H, V 
and U are endowed with the scalar products 

(m, v)h := {u, w)l2(o,r3), (m, v)v := {Lu, v)v'y, {u, -y)D(L) := {Lu, iu)L2(n,R3), 

respectively, and we denote by | • \h-, \ ■ \v and | • |d(l) the corresponding norms. 
Finally, for any integer fc > 0, we introduce the space Banach W*' of measurable 
vector functions u ~ (ui, M2, "3) defined in Q such that 

\u\y^k := ^ CSS s,up\d{d"u{t,x)\ < 00, (2.1) 
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where the sum is taken over < j < k and |q;| < 1. In the case k — 1, we will write W 
instead of W^. 

It is well known (e.g., see [20]) that problem (jl.ip . (|1.2|) is equivalent to the 
following evolutionary equation in H: 

ut + Lu + Bu^U{h + C), (2.2) 

where Bu : V ^ V is defined by Bu := B{u, u) with 

3 

{B{u, v), w)v' V — / Ui{diUj)wjdx. 

In the following, we will deal also with linear equations obtained from (|2.2p after 
replacing B by one of the operators M{u) and M*{u), where m S W is a fixed function, 
M{u)v = B{v,u) + B{u,v), and M*{u) stands for the formal adjoint of B(m) with 
respect to the scalar product on H: 




— UjdjVi)widx. 



Namely, let us consider the problem 

n+Lr + Vr = f, te/o = (0, 1), (2.3) 
r(0)=ro, (2.4) 

where B = M{u) or M*{u). 

Lemma 2.1. For any u eW, uq e H, and f e L^(/o, V), problem ([^ 
has a unique solution r G W{Iq,V,V') , which satisfies the inequality 

IHc(/o,H)+/ V?vAt+ j kt|y'dt<q|^|^o]|ro|l^ + |/li.(,„,y,), (2.5) 

J Iq J Iq 

where I stands for the closure of an interval / C K. Moreover, if f £ L^{Iq, H), then 
we have the inclusions y/tv G C{Iq, V), \ftv G L^{Io, U) and the estimate 

\^t^\liio,v)+ [ (V^Mt/)'dt<q|^|^„](|ro|l, + |/|i.(,„^^)). (2.6) 

J Iq 

Finally, if ro e V and f G L^{Io, H), then r G W{Iq, U, H) and 

Mc(7o,v)+/ IHD(L)di+/ |rt|l^dt<C[H^„](|roP,. + |/|i.(,„^H)). (2.7) 

J la J lo 

The proof of this lemma is based on a well-known argument, and we will not 
present it here. We refer the reader to the books [ini HO] for more general results 
on existence, uniqueness, and a priori estimates for solutions linear and nonlinear 
Navier-Stokes type problems. 
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2.2. Setting of the problem. Let us fix a function h £ i^(M_|_, H) and suppose 
that u e {R+ , 1/) n W solves tlie Navier-Stokes system 

ut + Lu + Bu ^ h, t > 0. 

Given a function uq Cz H and a sub-domain a; C 17, our goal is to find a finite- 
dimensional subspace £ C L^{ui,M.^) and a control ( G L^q^(R+, £") such that the 
solution of the problem 

ut+Lu + Bu = h + nC, ■u(O) = uq (2.8) 
is defined for all t > and converges exponentially to u, i.e., 

\u{t) - u{t)\H < C e-""^ fort>0, 

where C and k arc positive constants. 

Let us write L'^{Q,R^) as a direct sum L'^{Cl,R^) = H ® , where iJ^ denotes 
the orthogonal complement of H in L^. For each positive integer N, we now define 
A^-dimcnsional spaces E^^ C and Fjv C -ff as follows. Let {0^ | i S No} be an 
orthonormal basis in L^{ft, K'^) formed by the eigenfunctions of the Dirichlet Laplacian 
and let < /3i < /32 < • ■ • be the corresponding eigenvalues. Furthermore, let 
{et \ i S No} be the orthonormal basis in H formed by the eigenfunctions of the 
Stokes operator and let < ai < a2 < . . . be the corresponding eigenvalues. For 
each S No, we introduce the iV-dimensional subspaces 

En span{(/«, | i < N} C L^(r2,R^), Fn span{e, | i < N} C H 

and denote by Pn ■ i^(fi,M'^) — > En and Un ■ i^(fi,M'^) — > Fn the corresponding 
orthogonal projections. We will show that the required control space can be chosen 
in the form £m = X-^m, where x G C^(fl) is a given function not identically equal 
to zero, and the integer M is sufficiently large. 

Let us note that, seeking a solution of p.Sp in the form u — ii + wc obtain the 
following equivalent problem for v: 

vt+Lv + Bv + M{u)v^IlC, v{0)=vo, (2.9) 

where vq = uq^u{0). It is clear that it suffices to consider the problem of exponential 
stabilization to zero for solutions of (|2.9p . Thus, in what follows, we will study 
problem (|2J| . 

3. Main result for linearized system. We fix a function u G Lj^Qp(M+, V)r\W. 
In what follows, it will be convenient to write the control C entering (|2.9p in the form 
C = xPmi]: where rj takes its values in L^(f2,R'^) and x £ C^{il) is a nonzero function 
not identically equal to zero. Thus, we study the problem 

Vt+Lv + R{u)v = U{xPMV), (3.1) 
v{0) = vo, (3.2) 

where vq G H. We refer the reader to [TBI HO] for precise definitions of the concept of 
a solution for p. II) (and all other Navier-Stokes type PDE's). 
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3.1. Existence of a stabilizing control. Wc begin with the following result, 
which shows that one ean choose a finite-dimensional control exponentially stabilizing 
the zero solution for 

Theorem 3.1. For each vq G H and A > 0, there is an integer AI ~ C[\u\yj^\\ > 1 
and a control rf"^{vQ) £ L^(R+, Em) such that the solution v of system p.l[) . p.2p 
satisfies the inequality 

\v{r)?H < ^^i\vo?He-^\ i>0, (3.3) 

where Ki = C[\u\y^,^x] > is a constant not depending on vq- Moreover, the mapping 
Vq I— > r]^''^{vo) is linear and satisfies the inequality 

l^<'/^V^'(-o)L.(«^,^„)<-2K'ok, (3.4) 

for < A < A, where K2 ~ ^[\u\„ a (A-A)-i]- Fi'^^ol^Uj if '^o G ^> then 

Ht)\l < n:,\v^\le'^\ t>0, (3.5) 

where K3 = C^{^\y^, x] > does not depend on vq. 

To prove this theorem, we will need two auxiliary lemmas. For each r > 0, 
consider equation p.ip on the time interval Ir ~ (r, r + 1) and supplement it with 
the initial condition 

v{t) = wq. (3.6) 

Let us denote by Su,r(wo, rf) the operator that takes the pair (wq, '7) to the solution 
of (|3.ip . (|3.6p . By Lemma I^Tl the operator S^^r is continuous from H to C[Ir, H) fl 
L^{Ir, V) and from V to C(/^, V) n L^{Ir, U). We wih write Su^r{wo, v)it) for the 
value of the solution at time t. 

Lemma 3.2. For each N Cz N there is an integer Mi ~ C'[a,|m|w,-/v] ^ 1 such that, 
for every wq G H , one can find a control 77 G L'^{Ir, Emi) for which 

TlNSu^r{wo, vi){t + 1) = 0. 

Moreover, there is a constant depending only on \u\w {but not on N and r) such 
that 

< CMl. (3.7) 

Proof. Let us fix e > and consider the following minimization problem. 
Problem 3.3. Given M, N € N and Wq € H , find the minimum of the quadratic 
functional 

on the set of functions {v, rj) G W{Ir, V, V')xL'^{Qr, R^) that satisfy ((3Tl) and ((3^ . 

Theorem 15.21 implies that Problem 13.31 has a unique minimizer {v^^fje), which 
linearly depends on e H . Wc now derive some estimates for the norm of the 
optimal control ry*^. 
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To this end, the general theory of hnear-quadratic optimal control problems is 
applicable. We use here a version of the Karush-Kuhn-Tucker theorem (see Theo- 
rem [5TT]). Let us define the affine mapping 

F : W{Ir, V, V) X L2(g^,R3) H x L^{U, V'), 

{v, f]) ^ {v{0) ~ wa, Vt + Lv + M{u)v - U{xPmv)) 

and note that its derivative is surjective. Hence, by the Karush-Kuhn-Tucker theo- 
rem, there is a Lagrange multiplier {fi'^, q'^) e H x L'^{It, V) such that0 

J^(w% n^) - q") o F'{v^, f) = 0. 

It follows that, for aU (z, £) e W{Ir, V, V) x l?{Qr, R^), we have 

-(nwi;'(T -f 1), z[t + \))h + (z(t), ii')h + [ {zt + Lz + M{u)z, q')v',v dt = 0, 

(3.8) 
(3.9) 



2 ifj\ Ol- dt + (-n(xPA/0, Q')v'y dt = 0. 



Relation (|3.8p implies that q'^ is the solution of the problem 



qt - Lq" - B*(u)g' = 0, t G /r, (3.10) 
9"(r + l) = -2e-injvw'(T-|-l). (3.11) 



Furthermore, it follows from p.9p that 

277, = PMixql- (3.12) 
Combining ((3?T|) . ((3?T0)) . and ((3?T2]) . we derive 

-|((Z^ i;')H = (g^ w')ff + (g^ wDff 

^{q%n{xPMff'))H = \\PM{xctl\h. 

Integrating in time over the interval I^-, we obtain 

\PM{xq'{t))\\. dt = 2((g^(r + 1), v'{t+1))h - (9^(r), v^{t))h). 

Recalling now ((XTT|) . we see that 2(g'(r 1), w"(r -t- 1))^ = -e|g"(T + 1)||^ and 
therefore 

' \PM{xql\l^ dt + e\q\T + 1)1?, = -2(q^(r), «'(r))ff . (3.13) 



^The space H X L^(Ir, V) is regarded as the dual of H X L^(It, V), and the sign o stands for 
the composition of two Uncar operators. 
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Wc wish to use the truncated observabihty inequahty (|5.9p to estimate the right-hand 
side of (|3.13p . To this end, we take M — Mi, where Mi is the integer constructed in 
Proposition [231 Then, for every a > 0, we can write 

I \PMixq')\h + e\q^ir + < a\q^ir)\l + a-'\v^{r)\j, 

<^D^j^ \PMXq%^'it + a-^\v'{T)\]j. 

Setting a ~ (2_D^)^^, we obtain 

^ \PM{xq')\hAt + 2e\q\T + l)\l<ADM''H- (3-14) 

In particular, the family of functions {Puix't) I e > 0} is bounded in L^{Qr, M"^), and 
the family of solutions {u^ | e > 0} for problem p.ip . (|3.6p is bounded in L'^{It, V). 
It follows that the family {v^ \ e > 0} is bounded in L^{It, V). Thus, we can find a 
sequence e„ 0+ such that 

V'" ^IPMixq'") ^ V° in L\lr,EM), 
^ u° in L^iU, V), 
vt" v° in L'^ilr, V), 

where r/'^ £ L^ilr, Em) and € W{It-, V, V) are some functions. A standard limiting 
argument shows that is a solution of problem p.ip . p.6p with rj — if . Furthermore, 
it follows from (f3J4|) and ([3lT|) that 

|n^z;^(T + l)|2, = ^|g^(r + l)|2,<i^|u;o|2,^0 as e ^ 0. 

This convergence implies that IiNV^{T + 1) = 0. Furthermore, it follows from p.l4p 
that the function rf satisfies inequality p.7p with = ^D^. The proof of the lemma 
is complete. □ 

In view of Lemma 13.21 it makes sense to consider the following minimization 
problem. 

Problem 3.4. Given integers M,N > 1 and a function wq G H, find the 
minimum of the quadratic functional 

on the set of functions (w, i]) G W{Ir, V, V) x L^{It, Em) satisfying equations p.ip . 
p.6p and the condition IInv{t + 1) = 0. 

The following result shows that the control 77 constructed in Lemma 13.21 can be 
chosen to be a linear function of the initial state. 

Lemma 3.5. Let N > I be an arbitrary integer and let M be the integer con- 
structed in Lemma \3.2l Then for any wq € H Problem \3.4\ has a unique minimizer 
(i;"'^, 77"'^) e W{Lr, V, V) X Em)- Moreover, the mapping ^ (w"'^, r?"'^) 

is linear and continuous in the corresponding spaces, and there is a constant de- 
pending only on \u\w {but not on N and r) such that 

\r'^\hii^^E,,)<CMl- (3.15) 
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Proof. Let us fix e N, set 

WNilr, V, V) {v € W{Ir, V, V) I nNvir + 1) = 0}, 

and define X as the space of functions {v, rf) <S Wm{Iti V, V) x L^{It, Em) that satisfy 
equation p.ip . In view of Lemma [3.21 and the linearity of (|3.ip . X is a. nontrivial 
Banach space, and the operator A : X H taking {v,ri) to v{0) is surjective. Thus, 
by Theorem 15.21 Problem 13.41 has a unique minimizer (z;"''^, 77"'^), which linearly 
depends on Wq- Inequality (|3.15p follows immediately from (|3.7p . because the norm 
of r7"''^(wo) in the space L'^{Ir, Em) is necessarily smaller than the norm of the control 
function rj constructed in Lemma 13.21 □ 

Proof of Theorem \3.1[ Let us fix a sufficiently large > 1 and denote by Mi 
the integer Af constructed in Lemma 13.21 The main idea of the proof was outlined 
in the introduction: we use the operator 77"'"^ constructed in Lemma 13.51 to define an 
exponentially stabilizing control ij'^'^ consecutively on the intervals /„ = {n,n + 1), 
n > 0. Namely, let us fix an initial function vq € H and set[l 

r7"'^(f) = r7"'0(wo)(t) for t e lo- 

Assuming that 77"''^ is constructed on the interval (0,n) and denoting by v{t) the 
corresponding solution on [0, n], we define 

rj^^\t)=t'"iv{n)){t) for i £ /„. 

By construction, t]"'^ is an E'mi -valued function square integrable on every bounded 
interval. Moreover, the linearity of fj^''^ implies that 77"'^ linearly depends on vq- We 
claim that, if A^ S N is sufficiently large, then the solution v of system p.l[) . p.2p 
with rj = 77"'^ satisfies inequalities (|3.3p and p.5p . 
Indeed, it follows from (|2.6p that 

|z,(l)|^ = |5,,oK, t-'{vo))m <Ci\uM {\vo\l + \x\U^n)\t'°ivo)\hii,,E,.,)) , 
where we set E = Ei\j-^ to simplify the notation. Since Ili\jvi^) = 0, we obtain 

aN\v{l)\j, < Hl)\l < CiiuM {\vof + \x\Uin)\f''M\lm,^E,,)) ■ 
Using the continuity of (see Lemma [3T5|) and setting C'^ :~ ^xlxlioojfj), we derive 

Hl)\j, < «^Xl)Py < + C!,)\vo\^. 

Taking A^ so large that > e^{C[\u\y^] + C^), we obtain 

Hl)\l<e-'Ml. 

We may repeat the above argument on every the interval /„ and conclude that 

\vin + l)\j,<e-Mn)\'H- 



■^Recall that the operator r;"'"^ depends on TV. 
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By induction, we see that the solution of problem (|3.ip . p.2p with ij — rf'^ satisfies 
the incquahty 

Hn)\l<e-^-\v^\]j. (3.16) 
On the other hand, in view of (|2.6p . we have 

l^;Pc(/„,H) <qi«iw](kHll/ + lxli=^l^7"'"(«H)li2(/„,B,,)) < (qin|w] + c;)K'WlH- 

Combining this with (|2.5p . we see that v satisfies inequahty p.3p . 

We now prove (|3.4p . It fohows from (|3.15p and p.l6p that, for any A < A, we 
have 

< C^e^ E e(^-^)"|i;ollf < ^2\v^?h- (3.17) 

neN 

It remains to prove inequahty p.Sp . In view of (|2.6p . we have 

l\/*^«lc(/„,y) < ^[l«lw] (bHIff + 3|xli=c(o)ir"(«N)li2(,„,^„)) 

< (C[|fiU]+3C;)|t;(n)|2,. 

Combining this with p.l6p . we see that 

Hn+l)\l <CMr^)\l<Cie-^''\v^\l. 
In view of ()2.7p . for n > 1 we derive 

<^[|«lw]l«Nlv+3|xli==(o)l'?''"(«W)li^(/„,BM)<^2e-^^^^^ 
whence it follows that 

Ht)\l<C^e-^'\vo\l for t>l. 

Using again inequality p.7p . we conclude that p.Sp holds. The proof of the theorem 
is complete. □ 

3.2. Feedback control. In this section, we show that the exponentially stabi- 
lizing control constructed in Theorem 13. II can be chosen in a feedback form. Namely, 
let us fix a nonzero function x G C'o°(^) and denote by £m the vector space spanned 
by the functions x^ji j = 1, • • • , Af. Note that, due to elliptic regularity, the eigen- 
functions tjij ^'I'e infinitely smooth in £7, whence it follows, in particular, that £m is 
contained in (a;,R'^) for any cj D suppx- We will prove the following theorem. 

Theorem 3.6. Given ii e W and A > 0, let M = C^^^lw-M ^ ^ be the integer 
constructed in Theorem \3.1[ Then there are a family of continuous operators (t) : 
H — ^ £m and a constant k = C^^^^^x^ such that the following properties hold. 

(i) The function t i— > Kf^{t) is continuous in the weak operator topology, and its 
operator norm is bounded by k. 
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(ii) For any s > and Vq G H , the solution of the problem 

Vt + Lv + M{u)v = nK^{t)v, (3.18) 
v{s) = 1-0 (3.19) 

exists on the half-line Rs and satisfies the inequality 

e'^'-'^Ht)\l + j\^^'-'\Hr)\l + \vt{T)\l,)dT < Klvofn, t > s, (3.20) 
Moreover, if vq S V , then 

e^(*-^)b(t)Py+ f e''-^-^\\vir)\l^^^ + \vt{r)\l)dT < n\vo\l, t > s. (3.21) 

J s 

To prove this theorem, we wiU need two auxihary lemmas. Let us consider the 
foUowing problem. 

Problem 3.7. Given s > 0, A > 0, A/ e N and wq G H, find the minimum of 
the functional 

M^iv,v)--^ f e'\\v\l. + Hh)dt 

on the set of functions (v,ri) <E Wxi^s, V, V) x L\{Rs, i'^(il,K'^)) that satisfy equa- 
tion (13.11) and the initial condition 



v{s) = wo. (3.22) 

The following lemma establishes the existence of an optimal solution and gives a 
formula for the optimal cost. 

Lemma 3.8. For anyu e W and X > there is an integer M = > 1 such 

that Problem \3. 7| has a unique minimizer (v*,ri*). Moreover, there is a continuous 
operator Q^'^ : H ^ H such that 

M^{v:,7j:)^{Qfwo,Wo), (3.23) 
\Qf\ciH) < Ce^% (3.24) 

where C{H) stands for the space of continuous linear operators in H with the natural 
norm and C = C^\u\w->A > is a constant. Finally, Q'^ continuously depends on s 
in the weak operator topology. 

Proof Let X be the space of functions {v,r]) € Wx{Rs, V, V) x L\{M.s, L'^) that 
satisfy (|3.ip and endow it with the norm M^{v, ry)^^^. It is straightforward to see 
that A" is a Hilbert space. Using Theorem 13.11 with a constant A > A and the initial 
point moved to s, one can construct an integer M ~ C[|ii|^^ a] — -'^ such that, for any 
Wq G H and an appropriate control rj G -L|(Rs, -^a/), we have 

1^^)11, < ^ie-''^'-^^\wo\%, h(OlL < «2e-^(*-^)|^^o|lf, 
where v stands for the solution of p.ip . p.22p . Furthermore, by Lemma [531 '^e have 

e^'\v\'ydt < e^^^+'^ \v\ldt<C[\.a\„,x]e^^\v{r)\jj for any r > 0. 
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Combining the above three inequalities, we conclude that 

M^iv,ii)^ I e^*(IHy + HiOdt<%(i,_;,)-.,l^l^]e^1^^oll^. (3.25) 

J Ms 

It follows that X is nonempty, and the mapping A : X ^ H taking (w, 77) to v(0) 
is surjcctive. Thus, by Theorem 15.21 Problem 13.71 has a unique minimizcr {v*,rj*) = 
{v* (wq) , rj* (wq)) , which linearly depends on wq- 

We now prove p.23p and (|3.24p . It follows from p.25p that the mapping 

is a continuous bilinear form on H which is bounded by C2e^* on the unit ball. 
Therefore, the optimal cost can be written as p.23p . where Q^''^ is a bounded self- 
adjoint operator in H whose norm satisfies p.24p . 

It remains to establish the continuity of in the weak operator topology. To 
this end, it suffices to prove that 

{Q'^^^w, w) — > w) as s — > So for any w E H. (3.26) 

We will prove this convergence for sq > 0; in the case sq — 0, the proof is simpler. 

Let us fix T < So and denote by {v*,r]*) the optimal solution of Problem 13.71 
with s = T and wq ~ w. Note that, for any bounded interval / C Mr, the norms 
of the functions v* and r]* in the spaces W{I,V,V') and L'^{I,L^), respectively, are 
bounded uniformly in r. It is straightforward to see that the restriction of {v*,r]*) to 
the half-line Rs with s > r is the optimal solution of Problem 13.71 with wq = v*{s); 
cf. Lemma [3.101 below. Therefore, abbreviating = Q^' , we can write 

/CO 
e^\\vm\v + \vm\h)dt. 

Setting A^^{w) — \{Q'^v*{s),v*{s)) — {Q'^w,'w)\, for s > r we have 

\{Q'w,w) - {Q'°w,w)\ < A^(u;) + A;°(«;) + 

The third term on the right-hand side of this inequality goes to zero as s sg- 
Therefore convergence p.26p will be established if we prove that A^(u') -> as 
T, s — > So- To this end, note that the above-mentioned boundedness of v* implies that 

\v*{s) — 'w\l2 — > as T, s — > So. 

Combining this with the fact that the norm of Q'^ is bounded on finite intervals 
(see p.24p ). we arrive at the required assertion. □ 

Wc now consider another minimization problem closely related to Problem 13.71 
with s ~ 0. 

Problem 3.9. Given A > and vq G H , find the minimum of the functional 
N^iv, ry) / e^\\v\l + \v\h) dt + iQfv{s), v{s)) 



re^\K{t)\v + \vmh)dt. 

J Sn 
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on the set of functions {v, rj) e W{{0, s), V, V) x ^^((O, s), L'^{fl,M.^)) that sat- 
isfy p. II) . p.2[) . where M is the integer constructed in Lemma \S.8\ 

Theorem 15.21 implies that Problem 13.91 has a unique minimizer (v*, rj*), which 
is a linear function of vq € H. The following lemma is the dynamic programming 
principle for Problem 13.71 with s ~ 0. 

Lemma 3.10. Under the hypotheses of Lemma \3.8[ the restriction of {vQ^rj^) to 
the interval (0, s) coincides with (w*, 77*) and the restriction of (wg,77o) to the half- 
lineWs coincides with , ^7s)('^-'o(■s))■ 
Proo/. We will confine ourselves to the proof of the first assertion, because the 
second one is obvious. Let us define the function 



(4, W) 

Then we have 



{vl 7j:){vom for ie(0, s), 

iv*s,v:){vtis)){t) for teRs. 



Mo\z*,0*)=N^ivlr,:). 
On the other hand, the definition of (w*, rj*) implies that 

N^ivl r?:) < N^{iv*, 77o*)|(o,,)) < M,\v*, r^*), 
whence it follows that 

Mo^z*, 0;) < Mq^K, 

The uniqueness of minimizer for Problem 13.71 with s = implies that (zq, 9q) = 
(wq, tjq), and the required assertion follows. □ 

Proof of Theorem \3.6\ Step 1. It is straightforward to see that Problem 13.91 
satisfies the hypotheses of the Karush-Kuhn-Tucker Theorem 15. 1[ in which 

X = VK((0, s),Vy) X L^{{Q, s), L^{n,R^)), y = Hx L^{{0, s), V), 

J = N^, and : A" 3^ is the affine operator taking {v, if) to (v(0) — WQ,vt + Lv -\- 
M{u)v — W{xPmT])) . Hence, there is a Lagrange multiplier (/Ltg, qs) G HxL^{{0, s), V) 
such that 

(7V,^)'(^;,^^,:)-(M„g.)oF'(^;:,^y:)==0. 
It follows that, for aU z e W{{0, s), V, V) and ^ e L^{{0, s), L'^), we have 

e^\vl z)v At + 2{Qfv:{s), z{s))h + (z(0), 

(zt + Lz + B(u)2, gs)y',ydt = 0, (3.27) 

2/ e^'{Tj:,OL-dt+ f {-U{xPMO,qs)v',vdt^Q. (3.28) 
JO Jo 

In particular, we conclude from p.27|) that 

iq,)t - Lq, - M*{u)q, = 2e^*Lt-:(t). (3.29) 
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Since qs,vl G L^{{0, s),V), we see that dtQa £ L^{{0,s),V'), and therefore Qs € 
W((0, s),l/, V), whence it follows that Qs G C([0, s], H). Using again (|32Z1), we 
derive 

(?,(5) = -2Qf (3.30) 
On the other hand, relation p.28p implies that 

Tl's = \e-^'PM{x<ls). (3.31) 

In particular, rj*{t) is a continuous function of t with range in E^. Combining p.3ip 
and p.30p . we derive 

Recalling Lemma l3.10l and using the fact that s is arbitrary, we conclude that rj^ is a 
continuous function of time with range in Em and that 

rio{t) = -e-^'PM{xQ'^vl{t)) foraU t > 0. 

Thus, the optimal trajectory Uq for Problem 13.71 with s = satisfies p.lSp . where 

K^{t) := ~e-^\PMxQf- 

It is clear that K^{t) is a linear continuous operator from H to £m- Moreover, 
it continuously depends on t in the weak operator topology, because so does the 
family Q^^^ ■ Finally, it follows from p.24p that the norm of K^{t) is bounded by 
a constant depending only on A and \u\w- We have thus constructed a feedback 
control K^(t) possessing the properties mentioned in (i). Moreover, repeating the 
above arguments for Problem 13.71 with an arbitrary s > 0, we conclude that its 
optimal solution {v*,ri*) satisfies the relation ri*{t) = —e^'^*'PM{xQ\'^v*{t)) for t> s. 
Hence, if v{t) is the solution of problem p.lSp . p.l9|) . then 

{v{t)~e-^'PM{xQTv{t))) = «(t),7y:(t)) for t > s. 
Combining this with (|3.24p . we conclude that 

{Qfvo,vo)^ f {e^Mt)\'l^+e-^'\PM{xQ*~Mt))\l.)dt<Ce^^\vo\j,. (3.32) 

Step 2. We now prove inequalities p.20p and p.2ip for solutions of problem p.lSp . 
(f3TT9|) . Let us &x vq e H and denote by v the solution of (|3TT8l) . ([3?T9| . It is 

straightforward to see that the function z{t) = e^^^'^'>*v{t) satisfies the equation 

zt + Lz + M{u)z = ^z + K^(t)z. (3.33) 

Taking the scalar product of (j3.33p with 2z and using the uniform boundedness of 
the family K^{t), we derive 

^\z{t)\l + 2|z(t)|2, = X\z{t)\% + 2{K^{t)z, z) - 2{V>{u)z{t), z{t))H < C[\u\y,,x] \z{t)\]j 
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Integrating this inequality over the interval (s, t) with t > s, recalling the definition 
of z, and using Lemmas 13. 8[ 13.101 and inequality (|3.32p . we obtain 

e^'\v{t)\l + 2 / e^^\v{T)\^y dr < e^'^volj, + Ci [ e^^\v(T% dr 

J{s,t) •J{s,t) 

<e^^\vo\l + C, f e'^\viT)\ldT 

< C2e^'\vo\ji. (3.34) 

Furthermore, it follows from p.lSp that 

e^'\vtit)\l,<Cse^'Ht)\v- 
Combining this with (|3.34p . we arrive at the inequality 

e'Mt)\j, + 2 / e'^{\v{r)\l. + |z;,(r)|2.) dr < C,e'^\vo\jj, 

which is equivalent to (|3.20p . 

We now assume that vq e V. Taking the scalar product of (j3.33p with 2Lz, and 
using the Schwarz inequality and the uniform boundedness of the family K^{t), we 
derive 

^^\m\v + '2\m\l(L) - H^,Lz)h + 2{Kl{t)z, Lz) - 2{M{u)z{t), Lz{t))H 
< \Lz\l + C[\t.\„^x]\z{t)\l. 
Integrating this inequality over the interval (s, t) and using (|3.20p . we obtain 

J(s,t) J(s,t) 

<C^e^'\vo\l- (3.35) 

Furthermore, relation (|3.18l) implies that 

e'''\vml<C,e^Mt)\l(Ly 
Combining this with p.35p . we arrive at p.2ip . □ 

We conclude this section with a few remarks. 

Remarks 3.11. (a) Once a feedback control is constructed, it is easy to find a 
time-dependent Lyapunov function for the problem in question. Indeed, let U{s,t) be 
the operator taking vq H to v{t), where v stands for the solution of p.lSp . p.l9p . 
It is straightforward to check that the functional 

/oo 
\U{t,T)w\l2^a,M^)dT 



decays along the trajectories of (|3.18p . It is difficult, however, to write down this 
functional in a more explicit form. 
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(b) The operator Q*^''^ defining the optimal cost satisfies the following Ricatti 
equation: 

Q-{Q L(u) + L* {u)Q) - e'^'Q{WxPMX^)Q = -e^'L, (3.36) 

where L(u) = L + ]B(m). Since this equation does not play any role in this paper, we 
confine ourselves to its formal derivation. Let v be the solution of p.lSp , (|3.19p with 
vq G H and let T]{t) = —e~^*-PMxQ^v{t), where we set Q* = Q*^'^. By the dynamic 
programming principle (cf. Lemma l3.10p . the restriction of (w,?7) to the half-line Mi- 
is the optimal solution of Problem 13 . 71 with s = t and wq — v{t). Therefore, we have 
(cf. the equality in p.32p ) 

{Q-v{T),viT)) = [ {e^^v{t)\^y+e-^'\PM{xQ'v{t))\l.)dt. 

Differentiating this relation with respect to r and carrying out some simple transfor- 
mations, we obtain 

((Q^ - {Q^Uu) + L*iu)Q^) - e-^^Q^{UxPMX^)Q^ + e^^ L)v{t),v{t)) = 0. 

Setting r = s and recalling that v{s) — vq is arbitrary, we conclude that satis- 
fies (13311). 

(c) In the case of a stationary reference solution w, it is possible to give a rather 
sharp description of the dimension M for the feedback controller whose range depends 
on ii; e.g., see [4j [U [19]. In our situation, the range of the controller depends only 
on the norm of u, and its space dimension is determined by the integer Afi in the 
truncated observability inequality (|5.9p with a sufficiently large integer N. However, 
the feedback operator depends on time, and its image may be infinite-dimensional in 
time. It would be interesting to find out if it is possible to reduce the space dimension 
of the controller in our situation using further information about u. 

4. Stabilization of the nonlinear problem. 

4.1. Main result. Let us consider the nonlinear problem 

Vt + Lv + Bv + M{u)v^K^{t)v, t £ M+; (4.1) 
«(0) = vo, (4.2) 

where the operator K^{t) is constructed in Theorem l3.6l Given a constant A > 0, we 
denote by the space of functions z € C(K+, V) D Lf^^(M.+ , U) such that 

\z\z. := sup fe^*|z(t)Pv- + / e^l^MI^(i) dr^ < cx). 

t>o y j{t,t+i) J 

The following theorem is the main result of this paper. 

Theorem 4.1. Given u G W and A > 0, let M — C[\u\y^^x\ be the integer 
constructed in Theorem \3.b\ Then there are positive constants and e depending only 
on \u\w o-nd A such that for \vq\v < e the solution v of system (|4.ip . (|4.2p is well 
defined for all t > and satisfies the inequality 

b(i)ly <^e"^>olv fort>0. (4.3) 
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Proof. We will use the contraction mapping principle. We fix a constant i? > 
and a function vq e V and introduce the following subset of Z^: 

:= {zeZ^l z(0) = vo, < ^\vo\lr}. 

We define a mapping E : Z^ ^ C{R+, V) n ifo^(R+, U) that takes a function a e Z^ 
to the solution of the problem 

bt + Lh + M{u)b = K^b- Ba, t e M+, (4.4) 
6(0) - wo- (4.5) 

Suppose we have shown the following proposition. 

Proposition 4.2. Under the hypotheses of Theorem \^.l\ there exists > such 
that the following property holds: for any 7 £ (0, 1) one can find a constant e = > 
such that for any tjq G 1/ with |wo|y < e the mapping S takes the set into itself 
and satisfies the inequality 

|S(ai) - S(a2)|2:A < 7|ai - a2|2:A for all ai, 02 G Z^. (4.6) 

Thus, if \vq\v is sufficiently small, then the contraction mapping principle im- 
plies that there is a unique fixed point v G Z^ for S. It follows from the definition 
of S and Z^ that v is a solution of problem (|4.ip . (|4.2p and satisfies the required 
inequality (|4.3p . We claim that w is the unique solution of (|4.ip . (|4.2p in the space 
C(M+, F)nLfQ^(M+, [/). Indeed, if w is another solution, then the difference z = v — w 
vanishes at t = and satisfies the equation 

zt + Lz + B{z, v) + B{w, z) + M{u)z = K^{t)z. 

Taking the scalar product of this equation with z in H, carrying out some standard 
transformations (e.g., see [20]), and using the uniform boundedness of the feedback 
control K^{t) as an operator in H, we see that z = 0. Hence, to complete the proof 
of the theorem, it suffices to establish the above proposition. This is done in the next 
subsection. □ 

Remark 4.3. The hypotheses of Theorem [13] can be relaxed. Namely, it suffices 
to assume that the reference solution ii satisfies the condition 

sup(|m|l=c(q^) + |ut|i2(Q^)) < 00. 

T>0 

Indeed, as is proved in [7], the observability inequality (|5.8|) remains valid in this 
situation. It follows that the truncated observability inequality (|5.9|) . which is the 
key point of our approach, is also true. One can check that all the proofs can be 
carried out under the above weaker hypothesis. However, some calculations become 
cumbersome, and for the sake of clarity of the paper, we have imposed the more 
restrictive condition ii G W. 

4.2. Proof of Proposition l472l Step 1. We first derive an estimate for solutions 
of the equation 



zt + Lz + M{u)z = K^z + f(t), 



(4.7) 
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where / e Lf^^(M+,iJ). Namely, we will show that 

suJe''\zit%+ I e^1z(.)|2, d.) < Ci f |z(0)p^+sup / e^^^is^ds], 
t>o V J(t.t+i) J \ t>oj{t,t+i) J 

(4.8) 

where C\ = Cy^.^^^^ x^ is a constant. Indeed, recall that U{s,t) denotes the operator 
taking vq £ H to v{t), where v stands for the solution of p.lSp . p.l9p . By the 
Duhamel formula, we can write z as 

z{t) = C/(0, t)z{0) + [ U{s, t)f{s) ds. (4.9) 
J(o,t) 

Combining this with (|3.20p . we derive 

\z{t)\l = 2|t/(0, t)z{0)\j, + 2( f Pis, t)fis)\Hds^ 



2Ke-^^ziO)\l + 2Ke-^'(^J^ e'^^^^>\fis)\Hds^ . (4.10) 



'(o,t) 

Now note that, for any non-negative function c{t) and any A > 0, we have 

sup/ e(^/2).^(s)ds< / e^^/^^'c{s)ds^f] [ e^^/^^'c{s) ds 

t>oJ(0,t) J(0,+oo) ^^-^J(fc-l.fc) 

1/2 

\c{s)\'ds] 

/(fe-1. k) 



\J(k-l.k) 

oo . „ s 1/2 

<y^^-W2Kk-2)( e'^Ms)\'ds] 

k=l \j{k-l,k} 

/ r \ 1/2 

< C2 sup / e2^"|c(s)pds 

\t>Oj(t,t+l) 

Substituting this inequality with c{t) = |/(t)|ff into (|4.10p . we derive 

sup{e^'\z{t)\l)<2J\z{0)\j, + Clsup [ e'^^\f{s)\j, ds) (4.11) 
t>o V t>oJ(t,t+i) J 

On the other hand, it is easy to see that the analogue of Lemma 12.11 is true for 
equation (|4.7p . In particular, for any s > we have the estimates 



{t~s)\U{s,t)z^\\,<C^\\z^YH + ^Uirrndrj, (4.12) 
\U{s,t)zofy+ f |C/(s,T)2o|^(L)dT<C3f|zo|y+ / l/(r)|?^dT), (4.13) 

J{s.s+1) \ J(s.t) J 

where s < ^ < s + 1, and C3 > does not depend on s. Combining ()4.1ip with 
inequality (|4.12p in which zq ~ U{0, s)z{0) and < = s + 1, we obtain 

\z{s + < C3 f |C/(0, s)zml + f \f{r)\ldT) 

< C,e-'^ (\z{0)\j, + sup / e''^\f{T)\j,dT 

\ t>0 J(t,t+1) 
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Using now (j4.13p . for s > 1 wc derive 

W)\l+ ( |z(r)|^(^)dr<C5e-^^f|z(0)|l, + sup / e^^^/Wl^f dr) . 

J(s.s+1) \ t>OJ(t,t+l) / 

(4.14) 

On the other hand, it follows from (|4.13p that 

snp |z(s)|2,+ /■ \z{t)\I^^^At<cJ\z^\1+ [ |/(r)|l,dr). (4.15) 
o<s<i J(oa) V J(o,i) / 

The required inequality (|4.8p follows immediately from (|4.14p and (|4.15p . 

Step 2. We now prove that S maps the set into itself. Inequality (|4.8p with 
f{t) = —Ba{t) implies that 

|S(a)||. < Ci (\vofy + sup / e^^'\Ba{s)\j,ds) . (4.16) 

\ t>oJ{t,t+l) J 

Now note that l-Bajn < C6|a|v/|a|D(L), whence it follows that 

sup / e^^^\Ba{,)\\ds < C| sup / (e^->|^) {e^Ml^^^)ds < Cl\a\%,. 

t>Oj{t,t+l) t>oJ{t,t+l) 

Substituting this into (|4.16p . we see that if a € Z^, then 

|S(a)|2A < CrilvoW + \a\%.) <Cv{l + d\vo\v) \vo\v. (4.17) 

Setting -d = 2C7 and choosing e > so small that 6*7(1 + i9e) < we see that if 
l^oly ^ £: then S maps the set Z^ into itself. 

Step 3. It remains to prove that S satisfies inequality (|4.6p . Let us take two 
functions ai, 02 G and set a = ai — 02 and z = S(ai) — 2(02). Then the function z 
satisfies the initial condition z(0) = and equation (|4.7p with / = _Ba2 — -Boi. 
Therefore, by inequality (|4.8p . we have 

|S(ai) - S(a2)||. < sup / e'^'\Bai - Ba2\]jds. (4.18) 

t>0 J(t,t+1) 

Using a standard estimate for B{u,v) and the inequality |u||oo < C|M|v/|u|D(i), we 
derive 

\Ba1-Ba2\l = \B{aua)- B{a,a2)\ji 

< Cs{\ai\L'^\a\v + |a|L==|a2|y)^ 

< C9(|ai|v|ai|D(i)|a|y + |a| v|a|D(L) |a2| ■ 

It follows that 

/ e^^'\Bai^Ba2\lds<Cio{\ai\l, + \a2\l.)\a\l,. (4.19) 
J{t,t+i) 

Substituting (|4.19p into (|4.18p and recalling the definition of Z^, we obtain 
|S(ai) - S(a2)||. < 2i?Cio|wo|y|ai - a2||.. 

Choosing e > so small that < 7^, we see that if Iwolv < then (|4.6p holds. 

This completes the proof of the proposition. 
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5. Appendix. 

5.1. Karush—Kuhn— Tucker theorem. Let A" and 3^ be Banach spaces and let 
J : X M. and F : X ^ y he two continuously differentiable functions. Consider the 
following minimization problem with constraints: 

J{x) min, F{x) = 0. (5.1) 

We will say that x £ X is a local minimum for (|5.ip if F(x) — and there is a 
neighborhood U 3 x such that J{x) < J{x) for any x G U such that F{x) = 0. A 
proof of the following theorem can be found in |15| . 

Theorem 5.1. Let x £ X be a local minimum for ()5.ip and let the derivative 
F'(x) : X ^ y be a surjective operator. Then there is y* G y* such that 

J'ix)+y*oF'{x) = 0. (5.2) 



5.2. Quadratic functionals with linear constraint. Let X and y be normed 
vector spaces, let J{x, y) be a bounded symmetric bilinear form on X that is weakly 
continuous with respect to each of its arguments, and let A : X y he a continuous 
surjective linear operator. Given a vector y € y, consider the minimization problem 

J{x) — > min, Ax ~ y, (5-3) 

where J{x) ~ J{x, x). We will say that i G A" is a global minimum for (|5.3p if Ax = y 
and J{x) < J{x) for x & X such that Ax = y. The following result is rather standard 
in the optimal control theory, even though we were not able to find in the literature 
the statement we need. 

Theorem 5.2. Suppose that J{x) is non-negative and vanishes only for x ~ 0, 
and the set {x G A" : J{x) < c} is weakly compact for any c > 0. Then problem (|5.3p 
has a unigue global minimum x € X , and the function L : y ^ X taking y to x is 
linear. 

Proof. Existence. Let m > be the infimum of J on A^^{y) and let {xn} C 
A^^ly) he a sequence such that J(x„) — ?> to. Since the set {x £ X : J{x) < to + 1} is 
weakly compact, we can assume that {xn} converges weakly to a vector x £ X . Now 
note that 

< J{Xn — x) = J[Xn) — 2J[Xn,x) + J{x). 

Combining this with the weak continuity of J, we see that 

J{x) < liminf J(a;n) = rn. 

Thus, X is a, global minimum for (|5.3|) . 

Unigueness. Since the only point of X at which J vanishes is x = 0, a standard 
argument proves that J is strictly convex, that is, 

J{^H^) < 5('^(^i) + -^(^2)) forallxi,X2G A", 

and the equality holds if and only if xi ~ X2. This immediately implies that the 
global minimum is unique. 
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Linearity. Let y G y and z E ^^^(0). For all A > 0, wc have A{Ly ± \z) = y, 
and the definition of L implies that J{Ly) < J{Ly ± Xz). It follows that < XJ{z) ± 
2J{Ly, z) for all A > 0. Letting A go to 0, we see that 

J{Ly,z) = for all y€y,z€A-\0). (5.4) 

For a, b E y and a, /3 G M, let us set 

k := aLa + jSLb - L{aa + jSh). 

Then Ak = 0, and by (fO)) . we have J(fc) = J(fc, k) = aJ{La, k) + (3J{Lb, k) - 
J{L{aa + fib), k) ~ 0. It follows that fc = 0, and therefore L is linear. □ 

5.3. Truncated observability inequality. Wc first recall a well-known ob- 
servability inequality for the linearized Navier-Stokcs system. Let us fix a function 
u G L'^{It, V) n Wt, where Wr stands for the space of measurable vector- functions 
on Qr such that (cf. (f2TT|) ) 

|u|>v^ ;= ess sup \&ld^u{t,x)\ < oo, 

where the sum is taken over j = 0, 1 and \a\ < 1. Consider the problem 

qt- Lq-M*{u)q = 0, telr, (5.5) 
g(T + l)-gi, (5.6) 

where qi E H. By Theorem 2.2 in [M] (see also [7]), for any open subset w C there 
is a constant Cui such that 



\q{r)\ji<C^ I 191^2(^,83) di, 



(5.7) 



Since suppx n 7^ 0, the domain {a; e | |x(2;)| > p\ is nonempty for a 

sufficiently small p > 0. It follows from (|5.7p that 



^(t)!^^ < a, / kli.(.^,K3) di < c^^p-' / |x9li. dt. 

Thus, setting D'^ := C^^p^'^, for any solution of system (|5.5p . (|5.6p . we have the 
observability inequality 



\q{r)\l<D'^j^ 



xq{t)\l2dt. (5.8) 



The following proposition shows that if qi belongs to a finite-dimensional subspace 
of H , then the function xq on the right-hand side of (|5.8p can be replaced by Phiixq) 
with a sufficiently large Af . 

Proposition 5.3. For any N > 1 there is an integer Mi = (^[liiiw^,^] G No 
SMc/i that any solution q for system (|5.5p . (|5.6p wiit/i gi € i^Ar = UnH satisfies the 
inequality 



\q{r)\i,<D^ 1^ \PMAxqmhdt (5.9) 
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for a suitable constant depending only on x- 

To prove the proposition, we need the following lemma. 

Lemma 5.4. For any solution q of system (|5.5I) . (|5.6p with qi G Fj^, we have 

IX9WIhi(o,m^) dt<C f \xqit)\hin,R^) dt, (5.10) 

where the constant C depends only on N, f2, and \u\w- 

Proof We argue by contradiction. Suppose there is a sequence (g", u") £ F^ x 
{L'^{It, V) n Wt), with (|u"|w^) bounded, such that the solution (/" of the problem 

q,"-Lg"-B*(?i")g" = 0, t G (5.11) 

9"(r+l) = gl' (5.12) 

satisfies the inequality 

/ \xqlmdt>n [ \xq"\hdt. (5.13) 

Since the equations are linear, there is no loss of generality in assuming that = 1. 
The boundedness of (|u"|w^) implies that (9"zi") and (^"m") are bounded in L°°{Qr) 
for |a| < 1. It follows from Lemma [2Tl that the sequences (q") and (q") are bounded 
in L^(/r, D(i)) and L'^{It, H), respectively. Since the unit ball in a Hilbert space 
is weakly compact and the unit ball in L°"{Qt) is compact in the weak* topology, 
there is a subsequence of ((7",g",u") (for which we preserve the same notation), a 
unit vector qf" G Fn, and functions q°° G 1^(7^, D(L), iJ) and u°° G Wt such that 

q'^ qr in Fn, 

g" ^ 9°° hi L\lr,V), 
dtq'' dtq°° in L^I^, H), 

u!" u°° in L^{lr, H), 

did^ii'' in i°°(Qr), 

where j = 0, 1 and |a| < 1. Combining this with the boundedness of the se- 
quences (m") and (g") in the corresponding spaces, we can easily pass to the limit 
ill (|5.1ip . (|5.12p and derive the equations 

gt°°-Lg°°-B*(w°°)q°° ==0, t G I^, (5.14) 
q^{r + l)=qr. (5.15) 

Furthermore, since multiplication by x is a continuous operator in L^{Ir, FI^), we 
also have 

X9"^X9°° in L\It, H\n,R^)). (5.16) 
Therefore, passing to the limit in inequality (j5.13p as n — > c», we conclude that 

Xg°°|i.dt = 0. (5.17) 

Applying now the observability inequality (|5.8p to equation (|5.14p considered on the 
interval (t + r, r + 1) with < r < 1, we conclude that q°°{t) = for t < t < r + 1. 
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Since q°° S C(/^, V), we obtain q"^ — q°°{T + 1) = 0. This contradiets the fact that 
^oo ^ ^jj^ji- vector. The contradiction obtained proves that (|5.10p holds. □ 

Proof of Proposition 15.51 We use Lemma 15.41 to derive 

/ \xq\hdt< [ \PMixq)\hdt+ f \{l-PM)xq\hdt 

Jl-r Jlr JIt 

< I \PM{xq)\h dt + Pll I 1(1 - PM){xq)?m dt 

Jlr JIt 

< I \PM{xq)\l^dt + Pll I \xq\^I^dt 

JIt JIt 



< 



It JIt 

\PMixq)\hdt + pllC[N,\uU^] I \xq\l^dt. 



'It 

RecaU that /3j stands for the j*^ eigenvalue of the Dirichlet Laplacian. Choosing the 
integer M = Mi so large that /3MjC[Ar,|u|>v^] ^ "^(^ obtain 



1^ |xg|i.dt<2^ \PMAxq)\hdt. 



Combining this with (|5.8p . we arrive the required inequality (|5.9p . □ 
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